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Unitarity, a powerful constraint on new physics, has not always been properly 
accounted for in the context of hidden sectors. Feng, Rajaraman and Tu have sug- 
gested that large (pb to nb) multi-photon or multi-lepton sugnals could be generated 
at the LHC through the three-point functions of a conformally-invariant hidden sec- 
tor (an "unparticle" sector.) Because of the conformal invariance, the kinematic 
distributions are calculable. However, the cross-sections for many such processes 
grow rapidly with energy, and at some high scale, to preserve unitarity, conformal 
invariance must break down. Requiring that conformal invariance not be broken, 
and that no signals be already observed at the Tevatron, we obtain a strong uni- 
tarity bound on multi-photon events at the (10 TeV) LHC. For the model of Feng 
et al., even with extremely conservative assumptions, cross-sections must be below 
25 fb, and for operator dimension near 2, well below 1 fb. In more general models, 
four-photon signals could still reach cross-sections of a few pb, though bounds below 
200 fb are more typical. Our methods apply to a wide variety of other processes and 
settings. 
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I. INTRODUCTION 

The current era is dominated by hadron colliders, where small signals must be extracted 
from very large data sets. In order that new physics of an unfamiliar sort not be missed, 
it is important to consider a wide variety of possible signals that the experimenters might 
encounter. In this spirit, there has been considerable activity aimed at thinking broadly 
about reasonable non-minimal extensions of the standard model Higgs sector, of minimal 
supersymmetric models, and so forth. While there are strong motivations for each of these 
classes of models, the simplicity of their minimal versions is motivated mainly by aesthetic 
considersations. Moreover, the extra particles in non-minimal versons can lead to completely 
different phenomenological signals from those arising in the minimal versions. Given the 
baroque nature of the standard model, we would be unwise when addressing important 
issues in particle physics not to consider the possibility of particles and forces beyond the 
minimal set required. 

Considerable attention has been paid recently to hidden sectors that couple to the stan- 
dard model at or near the TeV scale. These include "hidden valleys" [H [21 E], new sectors 
with mass gaps and non-trivial dynamics, which lead to new light neutral particles, often 
produced in clusters and with a boost, and possibly with macroscopically long lifetimes. Hid- 
den valleys are especially natural hosts for dark matter, and indeed a class of hidden valley 
models |1] are a popular explanation for current anomalies in dark-matter experiments. 

Work on hidden sectors also includes a great deal of research on conformally invari- 
ant hidden sectors, dubbed "unparticles" in [5l EJ [7] (see also [HI |9]). New sectors with 
conformally-invariant physics (or at least scale-invariant physics, though there are no known 
examples of theories in four dimensions with scale invariance but without conformal invari- 
ance) can produce large missing-tranverse-momentum ("MET") signals, and can produce 
smaller, but potentially still dramatic, visible effects. However, the literature on this subject 
is full of contradictions, and many claims of interesting effects have been criticized. This 
has left the experimental community without clear guidance as to how to search for hidden 
sectors of this type. 

Our goal in this paper is to bring some clarity, through simple arguments, to a claim 
that large production rates for multi-particle final states can be generated through the three- 
point function of hidden sector operators that couple to the standard model. (Other work 
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emphasizing the importance of higher-point functions, often called "unparticle interactions," 
can be found in [TT[ |T2]. Additional subtle issues are addressed in [6l [3 [131 El US]-) We 
consider specifically the mechanism discussed by Feng, Rajaraman and Tu in [10], slightly 
generalized. In pUj it was pointed out that (for example) if a scalar primary operator O 
in the hidden sector couples to two gluons and also to two photons, and has a non-trivial 
three-point function {OOO), then the process gg — > 7777 can be generated. Because the 
form of a three-point function {O1O2OS) of primary scalar operators is precisely determined 
in conformal field theory in terms of the dimensions Aj of the three operators Oi, the 
kinematics of any process of this type is precisely known. (This is also true in some cases 
for three-point functions involving operators with non-zero spin.) In the case considered by 
[TU] . all kinematic distributions can be calculated in terms of the dimension and spin of O. 

Moreover, there is only one unknown parameter, the overall coefficient of the three-point 
function (equivalently the OPE coefficient connecting OO ^ O.) In [lOj it was pointed out 
that as of yet there is no known bound in four-dimensions on the size of this coefficient, 
and so it was suggested it could be arbitrarily large. Based on the limits from Fermilab on 
multi-photon events, it was claimed in [TU] that LHC production rates (at 14 TeV) were little 
constrained, and could range as large as 4 pb for ~ 1.1 and 8 nb (ten times larger than 
the ti cross-section) for Ao ~ 1.9. Given that four-photon backgrounds are tiny, and that 
the photons produced in this process have very high pt, this would be a truly spectacular 
signal by any measure. 

In this paper we throw some amount of cold water on this possibility. We first observe 
a simple-minded (and model-independent) unitarity constraint on any hidden sector, con- 
formal or not. Then we show how this specifically constrains conformally-invariant sectors, 
where explicit computations are possible due to the conformal invariance. After putting 
some experimental and theoretical limits on the size of the coupling between the two sec- 
tors, we apply this constraint specifically to the process pp 7777- For the specific case 
studied in [10], we find the maximum cross-section (for LHC at 10 TeV) is actually of order 
20 fb. When we generalize the scenario considered in [TU] by allowing the two gluons to 
couple to one operator Oi and the two photons to couple a different operator O2, we find 
that the maximum cross-section is anywhere from several pb, in the region Ai ~ 1.4 and 
A2 ~ 1.1, down to 30 fb or below for Ai + 2A2 > 5. 

Our methods can be applied more widely to various other processes. They will strongly 
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constrain four-lepton production through vector unparticles, for example, and any other 
similar process. 

As this paper neared completion some additional work on this subject appeared in p^[T7] . 
We believe that application of our methods would affect the conclusions of these papers. 
Also, in [17] production of multiple particles through exchange of two unparticles was con- 
sidered. While we do not address this issue in our current paper, there are additional and 
related unitarity bounds on this process which were not considered in [1^. It should also 
be noted that the authors of [T7] assumed in their calculation that there is no important 
four-point function among the hidden-sector operators, which is not universally true. 

The paper is organized as follows. We will explain our unitary bound in section II. After 
some general comments in section III about applications to unparticle sectors, we will show 
how to apply it to the specific case of gg — > 7777 in section IV. Section V will be devoted to 
obtaining a bound on the scale Ai characterizing the coupling between the two gluons and 
the unparticle sector. In section VI we will calculate the numerical bounds on pp — > 7777. 
We will comment on other possible processes in section VII, and state some conclusions in 
section VIII. 

II. A TRIVIAL UNITARITY BOUND 

We begin by pointing out an essentially trivial but rigorous unitarity bound that governs 
parton-parton cross-sections for hidden-sector production. The point, simply stated, is that 
no one process that involves the hidden sector can have a rate that exceeds the total rate for 
all such processes. 

This simple-minded and obvious point becomes useful when the total rate can be com- 
puted. Among the situations where this is possible is the case when the hidden sector is a 
conformal field theory to which the standard model (SM) couples via a local interaction. In 
this case the total cross-section is given by the square of a standard-model amplitude times 
the imaginary part of a two-point function of a local operator in the conformal field theory 
(recently given the name "unparticle propagator" [5].) Consequently, one may calculate the 
bound on the sum of all processes involving the hidden sector. 

Let us make a technically more precise statement of this unitarity bound. Suppose the 
interaction between the two sectors is governed by a local interaction, for example of the 
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form 

^^a^bO (1) 

where ipA,B are SM fields that create the SM partons A,B, and O is a gauge-invariant 
operator in the hidden sector that carries no SM charges. (We take O to be spinless for the 
moment, but our statements generahze for any spin.) 

We consider first a process AB —>■ X where X is a state in the hidden-sector Hilbert 
space. We will refer to the sum over all such states as AB {X}- Then the optical 
theorem assures that for center-of-mass momentum = q'^ + and center-of-mass energy 

a{AB ^ {X};s) = J](t(A5^X;s) 

X 

^ ImjAB ^ {X} ^ AB) ^ \{AB\^B^A\0)f lm[t{0\Oiq)Oi-q)\0)] 
s A'^^ s ' 

(2) 

Corrections to this last formula are smaller than the leading expression by a factor of order 
(s/A^)"^. We simplify notation by defining 

fAB ^ {AB\ijBi^A\0) ; Go{q;A)=z{0\O{q)O{-q)\0) (3) 

so that 

a{AB ^ {X};s) = ^^\f^^\Hm[Go{q; A)] (4) 

with s = q^. 

We are effectively assuming that the two sectors are weakly coupled to one another, so 
that the Hilbert space factors into a SM part and a hidden-sector part. This is true in the 
limit A — > cxD, and the corrections to this assumption should be small as long as momenta 
are small compared, naively, to 47rA. Actually, whether the condition involves 4'7rA or a 
somewhat smaller scale depends, as we will see, on the operator and on A, B. Also we have 
assumed here that any process generated by two separate couplings of the initial state to 
the hidden sector, such as considered in [17], is subleading compared to the effect of a single 
such coupling. If this is not the case, self-consistency problems arise, which we will not 
address here. 

Importantly, as emphasized by our notation, the two-point function of O that appears 
here is the complete two-point function, which includes all effects that depend on A from 



the interaction ([T]), along with any other interactions between the SM and hidden sectors. 
Let us define the two-point function of O in the hmit A ^ oo to be 

G^S\q) ^ \im Go{q; A) (5) 
The difference between this function and the full two-point function includes terms such as 

Go{q; A) = G'S\q) + ^G'S\qy^ J '^j'^ ^ 0\Mk)Mq - k)M-k)Mk - q)\0 > + . . . 

(6) 

as shown in Fig. [T| This particular type of correction sums as usual into a geometric series 



GS\q) 
l-Gg)(g)S(g) 



where 



s(g) = ^ / < ^\Mk)Mq - k)M-k)Mk - q)\0 > + . . . (8) 

as in Fig. [T} Other processes that connect the two sectors will also contribute to the full 
two-point function. 
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FIG. 1: The full two point function for O (filled line) differs from the conformal two-point function 
(unfilled line) by loops of standard model particles; these can be resummed as usual into a geometric 
series. 

Suppose we demand that the full two-point function Go{q', A) does not differ much from 
its A ^ oo limit Gq^ (q) — that is, that the interaction with the SM sector does not strongly 
alter the hidden sector in the energy regime of interest. (In particular, if the hidden sector 
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is conformal in the A — * oo limit, then we are demanding that it remain so to a good 
approximation.) Then any process such as AB P1P2 ■ ■ ■ + Xq, where Pi are SM particles 
and Xo is any hidden sector state, and where Pi are produced dominantly through SM- 
hidden sector interactions suppressed by 1/A to some power, can be bounded. In particular, 
this process will appear in the imaginary part of the full two-point function, suppressed by 
powers of 1/A to some power. The statement that the 1/A corrections to Go{q, A) are small, 
applied to its imaginary part, then implies that 

Y,<AB^P,P2---+X,) < J^|/^^|2im[Ga(g;A)-GgHg) 

< -r^\fAB\^lm[Go{q;A)] 



1 



- ^2S^Mlm [G'S\q)\ ^ a{AB ^ {X}) (9) 

The sum over Xq is over any subset of (and including possibly all) allowed hidden-sector 
states. The corrections to the last approximate equality vanish as A 00. Note the 
expressions in the first line are of higher order in 1/A than those in the last line, since by 
definition Go{q', A) —>■ Gq^ as A ^ 00. Therefore this is obviously true when A ^ q. But for 
LHC signals we will be interested in the consequences when q and A are not well separated. 

In English, the relations ^ state the following. The first inequality says that the process 
in question is found in the imaginary part of Go which does not appear in Gq\ since the 
latter contains only processes involving the hidden sector alone. The second inequality says 
that the difference between Go and Gq^ cannot be large, if conformal symmetry is valid. 
The third approximate equality restates that Go and Gq^ must be similar, so we may use 
either one. The final approximate equality comes from Eq. (|4]). The last two inequalities 
become equalities in the limit A ^ 00. 

It is crucial that the constraint ^ depends on q, or V^, the partonic collision energy, 
not directly on the collider energy -y/s. Thus, at a hadron collider, this constraint must he 
applied at all relevant values of 
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III. APPLICATION TO CONFORMAL HIDDEN SECTORS (UNPARTICLES) 

A. Conformal invariance must break down 

If the hidden sector is conformal, then (^^■'(g) is determined, up to a normahzation 
constant. The canonical normalization is taken so that in position space the time-ordered 
two-point function is l/(47r^x^)^ (up to contact terms at x = 0); any other normalization 
factor can be absorbed into A. The Fourier transform to momentum space yields 

w = - «)"-^ ■ (10) 

Our normalization is the same as that used in [5j , simplified by the use of Gamma-function 
identities. 

Suppose we want to use conformal invariance to predict something in the hidden sector. 
Then we must demand that any corrections to the two-point function are small compared 
to the two-point function itself, which then implies the bound (|9|. In particular, for any 
particular process (such as gg —>■ 7777, as we will consider below) in which only SM particles 
Pi are produced through the hidden sector, 

aiAB ^ {X} ^ P1P2 . . . Pn) < (t{AB ^ {X}) (11) 

In fact the bound is much stronger than this; the sum of cross-sections for all such processes, 
producing any standard model particles and hidden-sector states, is smaller than a{AB — *■ 
{X}). If conformal invariance predicts cross-sections that violate this condition, then it 
is conformal invariance itself that must be violated, and thus it cannot be used to make 
predictions. 

To illustrate the issues, let us consider a Lagrangian with three terms that couple the SM 
to the hidden sector through couplings to scalar hidden-sector operators, of the form 

5C = ^Oi^aV'B + ^(^2V^1^2 + ^(^3^3^4 (12) 

Here 5i = Ai -|- dimipA + dimipB — 4, and similarly for ^2, ^3. (For the moment we take all 
three operators Oi to be distinct; the case where the operators are related will be dealt with 
later. We also assume 5j > 0; we will discuss this assumption later. The standard model 
fields ipi, which create particles Pi, may or may not be different from one another; we make 
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no assumptions about them as yet.) Then, purely from dimensional analysis, we have 

.^A,^^xy,>).^{fy (13) 

where A^o is a constant calculable from conformal invariance alone and which depends only 
on Ai and on |/abP- Meanwhile, 

.(^B^P,;l)^^7Vp,(A,.A..A3)(g (14) 

Here, as emphasized by [TU], iVp. is a constant which is determined by the dimensions of 
the operators Oi. We will see we do not need its exact form. The OPE coefficient C123 for 
O1O2 —>■ O3 determines the normalization of the {O1O2O3) three-point function. Again, its 
value will not be needed for our discussion. 

These expressions are valid up to the scale s where conformal predictions break down. A 
sufficient condition for such a breakdown would be that a{AB Pf, s) ~ cr{AB {X}; s). 
If ^2 + (^3 > 0, as we are assuming at the moment, then a{AB —>■ Pf, s) grows faster with 
energy than a{AB — > {X};s). Thus there is always a scale Smax at which the expressions 



in Eqs. (13) and (14) become equal. At best, conformal invariance can be used to make 
predictions only up to this scale. At scales of order or larger than Smax there must be large 
corrections to the two-point function of Oi. When this happens, we can predict neither 
a{AB {X}) — which requires the two-point function directly — nor ai^AB Pi) 
— which is predicted using the special form of the three-point function, whose derivation 
requires that the two-point function of Oi be its conformal form. 

B. Motivation for studying gg — > 7777 

We must first decide what physical processes to study, which requires us to address some 



subtle points. The reader only interested in our results can jump to Sec. [IV 

We will focus on processes involving gauge bosons only. Our reasoning is the following. 
The largest effects from hidden sectors would come from low dimension operators. Scalar 
operators have the lowest possible dimensions, as is well known from unitarity bounds |18] . 
(See also jlHl 120] for other famous and important applications of these unitarity bounds.) We 



will discuss operators of non-zero spin in Sec. VII The only standard-model scalar operators 
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of low dimension are of the form (1) F^^F^'^ or F^^F^'^ for one of the standard model field 
strengths, (2) the Higgs boson bilinear WH, or (3) fHf, where / is a SM fermion doublet 
and /' is a SM fermion singlet. 

Large couplings of the form fHf'O break chiral flavor symmetries and are extremely 
dangerous, especially for the light quarks found in the proton. Without powerful symmetries 
or fine-tuning, these interactions will generically induce large and excluded flavor-changing 
neutral currents, through processes such as //' f'f, f /'7, etc., mediated via effects 
of the hidden sector. Conversely, suppressing flavor-changing neutral currents by choosing 
small couplings (i.e., choosing a very large value for A), reduces all cross-sections involving 
the hidden sector by factors of s/ to a positive power. We are skeptical that there exists 
an elegant model-building strategy that would permit operators to couple to the light quarks 
with A of order 1 TeV and A not far above 1 without risking large K-K mixing. Conversely, 
as A approaches 2, our bounds come into force. (Couplings of SM fermions to vector 
unparticles do not break chiral symmetries and are much more reasonable, but we are only 
considering scalar operators at the moment.) Consequently, it is far more natural that the 
initial state coupling should be to gluons. 

In the final state, fermionic couplings might have a role to play; for example, flavor- 
changing constraints on couplings to bottom and top quarks and to tau leptons are somewhat 
weaker, and one could imagine larger couplings of the heavier fermions to a hidden sector. 
We will discuss the possibility of a such flnal states in Sec. |VII[ 

Couplings to Higgs bosons are very interesting but are complicated by the relatively 
large mass of the Higgs and by its expectation value. Examples of these complications 
are described in [lH [15]. To avoid these complications in this paper, we assume that the 
couplings Who are not large, which in turn implies that the rates for producing Higgs 
bosons are small. In any case, Higgs bosons produced through a hidden-sector's three-point 
functions will lead mostly to multi-jet states, which have large backgrounds. 

For these reasons, in order to keep our presentation simple, we will focus on the process 
99 ~^ 1111- This case is nice both because it is conceptually straightforward, is a spectacular 
LHC signal, and was studied in some detail in [lOj. There are nevertheless some fine-tuning 
issues with the signal, which we discuss below. 
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C. A comment on the naturalness and fine-tuning 

On general grounds, when a theory has a low-dimension scalar operator O, fine-tuning 
is typically (but not automatically) necessary to avoid generating the operator O itself in 
the Lagrangian. This operator would then itself serve as a relevant perturbation of the 
conformal field theory and conformal invariance would be lost at very high scales. 

To avoid this, one would ask that any such operator transform under a global symmetry, 
so that its appearance in the Lagrangian is forbidden. For example, O might be a pseu- 
doscalar instead of a scalar, or it might transform with a minus sign under some other Z2 
transformation, or be part of a large multiplet under a continuous global symmetry, etc. 
However, these solutions are not entirely satisfactory since we must in general break this 



very symmetry to allow terms of the form (12). We might require that the standard model 



operator also transform under the global symmetry (for example if the Oi are pseudoscalars 
we can couple them to F^^^F'^'^, instead of F^yF^'^ as was done in [in]-) But this is not en- 
tirely satisfactory, because a three-point function among three scalar operators transforming 
under a Z2 symmetry must vanish, and more complicated symmetries which allow a three- 
point function cannot generally be realized among SM operators. For example, we cannot 
couple two gluons to an operator transforming under a Z3 symmetry without breaking that 
symmetry. 

We might also appeal to supersymmetry to prevent O from being generated with a large 
coefficient. In models where supersymmetry breaking in the hidden sector occurs at a scale 
which is low compared to the TeV scale, as can occur in models of gauge mediation where 
the hidden sector learns of supersymmetry breaking only through its coupling to the SM, 
supersymmetry can forbid the appearance of chiral operators in the superpotential, and thus 
restrict the operators that appear in the Lagrangian, down to a rather low scale. In this case 
conformal invariance would still be valid in the regime of interest. But this is not automatic 
and at the very least involves non-trivial model-building; see for example [2T] . 

Even if we solve the problem of generating O in the action, there is still the operator 
O'^O, which is usually a relevant operator for A significantly less than 2. (Note this operator 
as written may not itself be an operator of definite dimension, but it can be written as a 
linear combination of such operators, and one of them will generally have dimension less 
than 4.) The question of whether O'^O is relevant, and, if so, why it is not present with a 
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large coefficient, is analogous to the question of the small value of the Higgs boson mass. In 
order even to have a discussion about scalar operators with A well below 2, we must assume 
either that this coefficient is somehow unnaturally suppressed, or that it is protected by a 
very weakly broken supersymmetry in the hidden sector, as in [21]. (For interesting but not 
yet sufficiently powerful results regarding O'^O, especially where O has dimension less than 
2, see 

This particular problem does not arise for A > 2, where the square of the operator is 
generally irrelevant. (It has sometimes been erroneously suggested in the literature that 
scalar "unparticles" do not make sense for A > 2. But this is simply a misinterpretation of 
standard singularities which require standard operator renormalization. All conformal field 
theories contain such operators — for example, the square of the stress tensor.) Our results 
can be applied to such operators, but as we will see, the bounds that we obtain for such 
operators are on the verge of putting the signals out of reach of the LHC 

One may also ask about the coupling WHO, where H is the standard model Higgs 
boson. When the Higgs gets an expectation value, this inevitably would generate a breaking 
of conformal invariance [T^ [T3] . Again, if the conformal theory has an exact or weakly broken 
global symmetry that acts on this operator would be forbidden. (Meanwhile the operator 
W HO'^O is generally irrelevant.) In the models we consider below, any such symmetry is 
broken by the couplings to the standard model. But as long as the high-energy physics that 
generates these couplings does not directly couple the Higgs boson to the hidden sector, and 
a symmetry forbids WHO from arising well above the TeV scale, then any WHO term 
will be suppressed by an extra SM loop factor compared to the leading couplings between 
the two sectors, and will be sufficiently small not to undermine our assumptions. 

Thus to obtain gg —>■ '~f'~f'~f'~f from a conformally invariant sector requires quite a bit of 
work. But we will finesse all these issues, without further comment, in this paper. This is in 
order to address the specific phenomenological claims of [10], which assume implicitly that 
all these issue are resolved, but do not depend on the precise resolution. Also, although 
they are most easily explained in the case of scalar operators, our methods apply for any 
spin. At the end of this paper will briefly discuss more realistic settings, such as a three 
point function involving a vector operator V^, a scalar operator O, and its conjugate O^ . 
In this case the operator O could be a pseudoscalar, for instance, or carry some additional 
quantum numbers, and many of these problems would not arise. We emphasize, therefore. 
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that our results are very general and would apply with similar impact in many situations 
where there are no fine-tuning issues. 

D. A comment on the far infrared 

In general, conformal invariance in the hidden sector may not hold down to arbitrarily 
low energy. Indeed, we have just discussed various ways in which conformal invariance may 
be violated at low scales. Moreover, with the couplings that we consider, a truly conformal 
sector with very light particles can potentially induce new processes that have not been 
observed, or affect big-bang nucleosynthesis or other aspects of cosmology or astrophysics. 
For these reasons it may be that the hidden sector has a mass gap at some scale /x, which 
truncates all the branch cuts in Green functions of hidden-sector operators. (Examples of 
how this could occur appear in [121 [IH [IS]-) We will assume that any such fi is low enough 
that (1) it does not impact hidden-sector Green functions above a few tens of GeV, and (2) it 
does not cause any "hidden valley" signatures, where production of conformal excitations at 
high energy turns into hidden particles at the scale /i, which in turn decay to standard model 
particles on detector time scales, giving visible signatures [12] and completely changing the 
LHC phenomenology. We assume throughout this paper that any infrared effects do not affect 
the basic unparticle paradigm: that the hidden sector dynamics, for all observable purposes 
at the Tevatron and LHC, is conformally invariant and therefore predominantly invisible. 

IV. THE BOUND APPLIED TO FOUR-PHOTON EVENTS 

We now assume that the Lagrangian has couplings between the two sectors of the form 




where G° (a = 1,...,8) and F are SU{3) and [/(l)-electromagnetic field-strength ten- 
sors. For consistency, since the events we will study have energies far above the 100 GeV 
scale, we actually must couple the operator O2 to hypercharge bosons, with a coefficient 
{A2^ cos^ 6w)~^ ■ But for brevity we will ignore the associated 7Z and ZZ couplings for 
this paper. Although they contribute comparable three-photon and/or large MET signals, 
including them would not change the bounds that we obtain, which are in fact bounds on 
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the sum of the cross-sections for all these processes. Thus this omission is conservative, and 
simplifies our presentation. 

Note that we make explicit that Oi and O2 are distinct operators, potentially with 
Ai 7^ A2 and Ai 7^ A2. This need not be the case. They might be distinct operators with 
Ai = A2, or with equal Aj. Or we might take Oi = O2, as was assumed in [TU]; in this case 
we could assume Ai = A2, as in (TU], but we need not do so. In this sense our analysis is 
more general than that of [TO]. Indeed we will see the case they considered is much more 
strongly constrained than is the general situation. 

Now let us carry out our argument. Suppose, as we will obtain in the next section, that 
we have a lower bound on the scale Ai for given Ai. This is then an upper bound on the 
cross-section algg —>■ {X}; s) for producing anything in the hidden sector via the operator 
Oi. We could obtain from this a bound on the total hadronic cross-section a{pp —>■ {X}) 
by convolving this bound with the gluon distribution function in the proton. But this is not 
our goal. 

Instead, we turn to any particular process such as gg 7777, and require that it not be 
so large as to make preservation of conformal invariance impossible. In short, we require 

cr{99 1111; s) < a{gg {X}; s) {s < Smax) (16) 

But what V^max should we choose? 

To choose y/smax to be the collider energy would be too strong a condition. Most 
gg —* 7777 events at any collider will occur at energies far below the total collider en- 



ergy, and so y/smax need not be nearly so high. To determine the appropriate energy, we 
must compute the four-photon cross-section as a function of s, under the assumption of 
conformal invariance, and see where it is large. Then we should choose Smax so that the 
great majority of the 7777 events will be produced at energies below this value. 

For example, we might reasonably demand that a certain fraction ( of the gg 7777 



cross-section must occur below the scale y/smax- That is, we define Smax by 

d(r{gg im) ^ ^ r da{gg 7777) 
ds Jq ds 

where s is the square of the collider center-of-mass energy. To require C = 1, and therefore 
Smax = would bc far too strong, as noted above. If we instead took C = \ then we would 
effectively be demanding, typically, that the peak cross-section for gg 7777 occurs at 
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Smax, right where conformal invariance is breaking down. In this case, none of the predictions 
(cross-section or kinematic distributions) of [10] would be at all reliable. For this reason we 
view C = ^ as unreasonable. We therefore take C = | as a conservative choice. This should 
ensure that the prediction for the rate and differential distributions for gg 7777 are 
given to a rough approximation by conformally invariant calculations, and are not beset 
with model-dependent effects beyond roughly the 30%-50% level. 

Importantly, assuming only that conformal invariance has not been violated, we can 
determine Smax in a completely model-independent way that depends only on Ai and A2. 



From Eq. (14) (with Ai = 61 and A3 = A2 = §2 in the case at hand), we know the precise s 



dependence of the cross-section, up to constants that factor out of the condition in Eq. (17). 
Defining the gg luminosity function as usual by 



dr 



I dyfg{V^ey)fg{V^e-y) (18) 



(where t = s/s) and substituting from Eq. (14), we have, for ( 



I 



^^ ^^aai'^) ^Ai+2A2-i ^2 digger) 

dr 3 ./o dr ^ ' 



Notice that all dependence on C122, ^7777(Ai,A2) and Aj factors out of this expression. 
Thus our choice of s^axi once we have chosen a fixed C, depends only on Ai + 2A2, and 
largely scales with s (up to the slow variation of Lgg through the evolution of the gluon 
distribution function.) Table |l] shows \/smax for a 10 TeV LHC and various choices of 
Ai + 2A2. 

At this point we should mention that throughout this paper our numbers are produced 
using the (outdated) CTEQ5M pdfs [23j. This is purely for technical reasons of calcula- 
tional speed. Results obtained from more modern pdfs differ by significantly less than other 
systematic errors in our calculations. We have explicitly checked in several cases that our 
numbers do not change significantly with the MSTW08 pdf set |2l]. The errors on Smax 
from uncertainties in the gluon pdfs and the appropriate choice of factorization scale are 
estimated at approximately 5 percent. This is smaller than the dominant source of uncer- 
tainty, which arises from the choice of C that defines s^ax- We will have more to say about 
this uncertainty after we present our results. 

Now let us return to the process of obtaining a bound. The bound arises from the fact 
that oi^gg —>■ {X}; s) is precisely known, except for an overall constant normalization, which 
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Ai + 2A2 


3.0 


3.5 


4.0 


4.5 


5.0 


5.5 


6.0 


V Synax (in TeV) 


1.2 


1.7 


2.2 


2.7 


3.1 


3.4 


3.7 



TABLE I: Values of ^ 

Smaxi s 10 TeV LHC, for various choices of Ai + 2A2. 



depends only on Ai and is proportional to 1/A^^\ If Ai is bounded from below, Ai > A^*"', 
then a{gg {X}; s) is likewise bounded from above, at all s, by ai^gg {X}; s; A™"). 

To understand what this means intuitively, we have plotted CT{gg — > 7777; s; A^™) and 
'^{99 ~* {^}! -5) in figures |2] to [5| for several different choices of Ai and A2. The total 
hidden-sector cross-section (j{gg — > 7777; s; A^™) is normalized to saturate the bound on 
Ai that we will obtain later; however for the moment the shape matters more than the 
normalization. The normalization of the 7777 cross-section cri^gg — > {X}; s) is chosen so that 
it does not exceed the total hidden-sector cross-section at any below y/smax, whose value 
is indicated by a vertical line. Because of the rate with which the gg luminosity decreases, 
a{gg —>■ 7777; s) initially increases with energy, until the rapid decrease of the gg luminosity 
at high s overwhelms the rising partonic cross-section. Meanwhile, crlgg {X}; s) decreases 
rapidly everywhere. Because of this, we can see by eye that Smax must be taken quite large, 
typically of order 1-4 TeV. (This confirms that for gg — > 7777 we can neglect any effects 



from an infrared scale fx of the sort discussed in Sec. Ill D ) Also, we can see by eye that 
(T{gg 7777; s) is always vastly less than cr{gg {X};s), because of the shapes of the 
two curves, until s is very close to Smax- 

As we noted earlier in our more general discussion, dimensional analysis always assures 
that the ratio (r{gg 7777; s)/cr((7(7 {X}; s) grows with energy, as long as conformal 
invar iance is applicable. Therefore 

(^{99 7777; s) ^ (r{gg 7777; w) ^ (^{99 



a 



< 



7777; Sr 



(20) 



^ {X}; rs) a{gg ^ {X}; w) a{gg ^ {X}; s„ Ar") ' 
Unitarity requires the last expression be less than one, and writing this condition in terms 



of the constant coefficients appearing in the formulas (13) and (14) for the cross-sections, 
we obtain 

-2A2 



|Ci23rAr^777(Ai>A2)A2"'^' « (w)-^"^Xo(Ai 



(21) 



Notice all Ai dependence factors out of this bound. 
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dcr/dE 



Ai=].05 , A2=1.05 



O.OOlO 



0.0000 




0.001 




FIG. 2: For Ai = 1.05, A2 = 1.05, the differential cross-sections (in pb/GeV) versus energy E = 
(in TeV) for all production processes involving the hidden sector (thick curve) and for four-photon 
production (thin curve). The right hand plot is the same as the left hand plot, but on a log scale. 
The total hidden-sector cross-section is normalized by our bound on Ai, and the four-photon cross- 
section is normalized so that it satisfies unitarity, by not exceeding the total for any s < Smax- Our 
estimate of ^/smax-, determined as explained in the text, is indicated by the vertical line. 

Finally we may obtain a bound on the total cross-section for gg —>■ 7777, namely 

dLggir 



Jo 

« iVo(Ai) Af^Hsn^ax)-'^' s^^^'^^-' [\r 

Jo 



dr 



-T 



Ai+2A2-1 



'^-^m(''") ^Ai+2A2-l 



< 



iVo(Ai 



2A2 .1 




dr 



^-^gg(^) ^Ai+2A2-i 



(22) 



■ [A mm] 2 / \ £ / , , 

This is the formal expression of our main result. 

Notice that our bound only depends on the collider energy s, on the dimensions Ai and 



A2, on Smax/s (determined using Eq. (19) by Ai and A2,) on Nq (which is 

-sin(7rAi) r[2 - Ai] 



A^o(Ai) 



(47r)2Ai-2 



r[Ai 



(23) 



for a gg initial state,) on the known gg luminosity, and finally on A^*" (which we must 
separately determine using theoretical and experimental constraints.) All dependence on 
A2, N^^^^ and C122 has vanished. If we know A™" as a function of Ai and perhaps A2, we 
can obtain a bound that is model-independent and depends only on Ai and A2. 
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A,=1.35 , A2=1.05 




FIG. 4: Same as Fig. ^ but with Ai = 1.45, A2 = 1.45. 
V. OBTAINING BOUNDS ON Ai 

Our only remaining task is to determine A^*". Once we have it, we can compute the 
bound on the gg — > 7777 cross-section. 

We apply two main considerations for constraining Ai. The first is that if Ai is too low, 
then not only is the rate for the invisible process a{pp {X}) very large, the observable 
process a{pp j + {X}), where j is an initial-state jet, becomes comparable to the standard 
model rate for jet plus missing transverse momentum (MET). Contraints from Tevatron, 
mainly from the CDF study |2Sj, put strong contraints on Ai for low Ai. 

A second constraint on Ai comes from the fact that the coupling of O to gluons itself 
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dcrl dE 




FIG. 5: Same as Fig. |2l but with Ai = 1.85, A2 = 1.85. 

induces corrections to G^q . We must assume these are small if we are to use conformal 
invariance to make predictions regarding gg — > 7777. Either such predictions are impossible, 
invalidating the approach of [10], or Ai must be larger than some minimum. This puts 
moderate constraints, of order 1.5 TeV or larger, which are relevant for larger Ai where the 
experimental constraints are weakest. 



A. Bounds from Tevatron measurements of monojet events 

Given a known partonic cross-section for a hidden-sector process, it is straightforward to 
compute the rate for jets-plus-MET where the jet(s) only arise from the initial state. One 
might ask whether emission from the final state could possibly compete with, and perhaps 
interfere with, this process. The answer regarding interference is "no" ; once the hidden state 
has been produced, it is color neutral, and any final-state radiation must be color singlet, 
requiring at least two jets be emitted. Similarly, in the model we are considering, the largest 
interactions between the two sectors involve irrelevant couplings, so any final state radiation 
process is small at low energy, and is either too small to observe or would show up as a large 
tail at high energy. Since no such tail is observed at Fermilab, we assume any final-state 
radiation of jets cannot affect the limits which we will now obtain. 

For a conformal hidden sector produced through gg, the rate is entirely fixed by Ai and 
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Ai. For qg — > q{X}, we find, at leading order 
-{qg^q + {X}) = Cj2Bny<2Fi 



,-1 + Ai + n, -- + Ai +n, 



2 ' '2 ' s + 2\/ s-pT 

where p-r is the transverse momentum of the jet, 



(24) 



and 



(27rA)2Ai 3(3 + 2Ai)r[2 + 2Ai] p^^Vll - W/s) 



50 = (-3 - 2Ai + 12A2 + 8A^) (2 - 'ivrls) (26) 

51 = -2(Ai - 1) (3 + 8Ai + 4A?) (4 - 'ivT^js) (l - ^WT^) (27) 
^2 = 12(Ai - l)Ai(3 + 2Ai) (l - ywTs)' (28) 
^3 = -8(Ai - l)Ai(Ai + 1) (l - ^J^VTT^ (29) 



(The reader may compare our result with the literature, see for example [26], in the Ai ^ 1 
limit.) This is the dominant process at high energy at the Tevatron. There is also the 
process gg but this is smaller in the energy range of interest at the Tevatron and 

we neglect it. If we included it, our lower bounds on Ai would be stronger. 

The CDF experiment [23] has published results on monojet events, in the context of a 
search for extra dimensions, and a public webpage with additional information and plots 
is available [27j. Early results from DZero [28], with much lower statistics, have not been 
updated; we will not use them in our analysis. The CDF study uses two sets of cuts, a loose 
set for a model-independent search, and a tighter set optimized for an extra-dimensions 
search; we use the former. The data is available in plot form, though not in table form; we 
have extracted the data directly from the plots, introducing a moderate amount of systematic 
error in the process. Demanding that the process qg not be easily visible above 

the error bars of the plots in [25] puts a limit on Ai for any given Ai. 

Through this requirement we find limits on Ai shown in boldface in Table [Tlj There 
are substantial systematic error bars on our results. First, we have not included the K- 
factor from loop corrections, or the process gg 5'{-^}; doing so would give a slightly 
stronger bound. Second, we are not able to include experimental efficiencies and effects of 
jet energy scale uncertainties; doing so would give a slightly weaker bound. Furthermore, our 
computation is done at leading order, for which the jet transverse momentum p^^* and the 
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MET are equal. However, both additional jet radiation and jet mismeasurements contribute 
in the data, so these are not in fact equal, and thus when we extract a limit on Ai it is 
inherently ambiguous whether we should use the experimental distributions of da/d{MET) 
or d(7/d{pl^^) (and neither is accurate beyond leading order.) Crudely, we estimate that the 
errors on our determination of Ai are of order 10 percent, which turns out to be a subleading 
uncertainty compared to that stemming from the ambiguity in choosing Smax- 

As final comments, we note that for pl^* of this size, the cross section for da/dpl^^ involves 
an integral over that is insensitive to low g^. In other words, our limits on Ai are insensitive 
to any low-energy cutoff fi. Also, the reader may observe that our calculations do not suffer 
from the well-known singularity at A — 2 which indicates the need for renormalization. 
This is because our results depend only on the imaginary part of Go- AH of our results are 
smooth as A passes through 2. 



B. Bounds from preserving conformal invariance 

We noted earlier that in a conformal theory perturbed by an interaction of the form 
Eq. ([T]), there is an irreducible effect that causes Go{q] A) to differ from its conformal form 
Go\q), given by Eq. ^ and shown in Fig. 1 

At leading order, the QCD interactions of gluons play no role, and so we may treat them 
as a system of free massless particles — a conformal field theory. Thus our calculation is a 
specific example of a more general issue: if we have two conformal field theories / and J, and 
we couple them through an irrelevant operator OjOj with coupling 1/A^~'^^~^-', where Oj 
(Oj) is a scalar operator in conformal sector / (J), then this coupling leads formally to a bad 
breaking of conformal invariance at some high scale Mmax- More precisely, either conformal 
invariance is badly broken, or the pointlike coupling OjOj develops a non-pointlike structure 
due to new physics at some scale at or below M^ax- Either way, the approximation that 
one has two conformal field theories coupled by a pointlike operator must break down. 

What is an estimate for Mmax"^ With conventionally normalized operators Oj and Oj 
one might naively guess through naive dimensional analysis that M^ax ~ 47rA. With the 
normalization used in the unparticle literature (which sets the conventions for our definition 
of Aj in this paper), this is essentially correct. 

However, the standard model operator G'^^^G'^^^ is not a conventionally normalized 
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operator of dimension 4, because it contains derivatives. One may easily check that these 
produce additional factors of 27r (just as is expected in naive dimensional analysis) leading 
to a (27r)^ enhancement relative to the two point function of a conventionally-normalized 
operator of dimension 4. In addition, there is a factor of 8 = 3^ — 1 from the sum over 
colors. Altogether this means that, for the normalization of Ai given through the use of the 



action Eq. (15), which is the same as used by Feng et al. in |10j, the breakdown of conformal 
invariance occurs well below 47rAi. This is significant because in the literature one often 
sees discussion of taking Ai ~ 1 TeV, which may cause conformal invariance to break down 
within the range of energies accessible at LHC For our current problem, since the peak of 
the gg 7777 cross-section occurs at energies typically greater than 1 TeV (see Table [l] 
and Figs. |2]-|5|, this problem is severe. 

More precisely, the momentum-space two-point function of Gf^^G^" is quartically diver- 
gent, and there are underlying quadratic and logarithmic terms; renormalization removes 
these divergences but leaves their finite contribution ambiguous. However the imaginary 
part of the two-point function is unambiguous, arising from a finite In q term. When this 
imaginary part makes an order-one correction to G^^\q), conformal invariance is unambigu- 
ously breaking down. 

Even more precisely, we can see from Eq. ([t]) that we can no longer trust conformal 
invariance once \GQ\q)T,{q)\ is of order 1. As we have just noted S(g) is subject to renor- 
malization ambiguities, and for the same reason, so is Gg^(g) if Ao > 2. But the imaginary 
parts of S and Gq^ are not subject to such ambiguities. Noting 



Gg)(g)S(g)| > 



Im[Gg^(g)]lm[S 



we choose to apply an extremely conservative consistency condition, namely 



Im[Gg)(g)]lm[S(g)] 



< 1 



(30) 



(31) 



for any s < s^ax- This then gives a conservative lower bound on Ai. 
Explicitly, we find, in the notation of Eq. (|8|), 



a b 

Keeping only the finite imaginary parts, our consistency condition becomes 

,2 \ Al 



Im[Gg^(g)]lm[S 



sin(-7rAi)r[2 - Ai] 2 / g' 
^ (47r)2A-2r[Ai] ^ ^ 



< 1 



(32) 



(33) 



23 



for < Smax- Here the important prefactor of 8 counts the number of gluon states. This 
condition in turn imphes a lower bound on Ai. 

The uncertainties that arise here stem mainly from the ambiguity in the criterion chosen. 



For example, suppose we replaced 1 on the right-hand side of Eq. (33) with ^? This would 
only change Ai by (2)^/^^^, and strengthen our final bound by exactly a factor of |. This is, 
again, smaller than the uncertainty in our bound that arises from the ambiguity in defining 



As a final comment, we note that an analogous argument applies for many other standard 
model operators, including those with higher spin, putting similar lower bounds on the scale 
A. We are not aware of this constraint being accounted for elsewhere in the literature. 



A2 

Ai 


1.05 


1.15 


1.25 


1.35 


1.45 


1.55 


1.65 


1.75 


1.85 


1.95 


1.05 


9.19 


9.19 


9.19 


9.19 


9.19 


9.19 


9.19 


9.19 


9.19 


9.19 


1.15 


5.18 


5.18 


5.18 


5.18 


5.18 


5.18 


5.18 


5.18 


5.18 


5.18 


1.25 


3.19 


3.19 


3.19 


3.19 


3.19 


3.19 


3.19 


3.26 


3.43 


3.60 


1.35 


2.11 


2.11 


2.24 


2.43 


2.62 


2.80 


2.98 


3.15 


3.31 


3.47 


1.45 


1.76 


1.95 


2.13 


2.31 


2.48 


2.65 


2.81 


2.96 


3.10 


3.24 


1.55 


1.68 


1.85 


2.01 


2.17 


2.32 


2.47 


2.61 


2.74 


2.87 


3.00 


1.65 


1.59 


1.74 


1.89 


2.03 


2.16 


2.29 


2.41 


2.53 


2.65 


2.76 


1.75 


1.50 


1.64 


1.77 


1.89 


2.01 


2.12 


2.23 


2.34 


2.44 


2.54 


1.85 


1.42 


1.54 


1.65 


1.76 


1.87 


1.97 


2.07 


2.16 


2.25 


2.34 


1.95 


1.34 


1.45 


1.55 


1.65 


1.74 


1.83 


1.92 


2.00 


2.08 


2.16 



TABLE II: The minimum values of Ai (in TeV), as a function of Ai and A2, allowed by the 
experimental constraints from monojets and by the theoretical constraint that conformal invariance 
be preserved below Smax for the corresponding Ai, A2; see Table |I} Values constrained by monojet 
data are shown in boldface. 
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C. Summary of the bounds on Ai 

The bound we obtain from the more powerful of these two constraints, as a function 
of Ai and A2, is shown in Table [Tlj The constraint from jet-plus-MET measurements at 
the Tevatron is most powerful at small Ai, while the constraint of conformal invariance is 
the dominant effect at larger Ai. Notice that the conformal invariance constraints give a 
bound that becomes stronger as A2 increases, for fixed Ai. Note also that the bound never 
dips below 1 TeV. One should also keep in mind that bounds on monojets at Fermilab are 
probably stronger now than those which are currently published. The published CDF study 
[25] relies only on 1.1 pb~^. Though it is systematics-limited, it appears that some of these 
systematic uncertainties are data driven and will have decreased with higher statistics. 



VI. BOUNDS ON pp 7777 AT THE LHC 



With the bounds on Ai from Table pl|, we may now obtain bounds on a(pp 7777) using 
the condition from earlier sections. First we obtain bounds based on our central values and 
naive tree-level results; then we discuss their uncertainties. 



A. Bounds in the model of Feng, Rajaraman and Tu 

Let us consider first the particular case studied in (TU], where Oi = (92) = ^2 and 
Ai = A2. Because of the equal Aj, the processes gg gggg, gg 9911, and gg 1111 all 



have the same energy dependence, so unitarity constrains their sum, generalizing Eq. (22): 
a{pp gggg; s) + a{pp gg-fr, s) + 7777; s) 



^ '.[A— l^J VwJ Jo dr ^ 



(34) 



All processes listed here proceed through the hidden sector; QCD contributions to gg gggg 
are of course not to be included. 

To go further, we use the fact that the amplitudes for these processes are identical (since 
neither electromagnetic nor strong interactions enter the calculation at leading order); one 
may view the calculation as taking place in U (3) instead of SU (3)-color, with the photon 



25 



Ai = A2 


1.05 


1.15 


1.25 


1.35 


1.45 


1.55 


1.65 


1.75 


1.85 


1.95 


Max fj [maximal interference] (in fb) 


10.92 


19.26 


21.79 


15.63 


5.35 


1.98 


0.81 


0.34 


0.14 


0.07 



TABLE III: The maximum allowed values, in femtobarns, of the cross-section for pp — > 7777, as a 



function of Ai = A2, assuming Oi = O2 and Ai = A2, as in [10 . In this case — see Eq. (34| — 
both pp gggg and pp gg'J^ contribute to the unitarity bound. Since we have not performed 
the calculation directly we simply assume maximal interference among diagrams; the true bound 
obtained from such a calculation would be stronger, probably by a factor of 1.5-2. 

being the ninth gluon. The only non-trivial aspect is interference, which could be precisely 
computed, but we will only estimate. 

Label the gluons with an index a = 1, ... 8, with a = 9 for the photon. Label the matrix 
element for gg —>■ g'^g'^g^g^ as M.ab{ki, k2, k^, k^). Only the sums kij = ki + kj enter the 
amplitude. Then Mab = F{kl2, ^34) + F{kl^, k24)Sab + F{kl^, k^^jSab- Also for a = 6 there is 
a reduction in phase space by 3, due to Bose statistics. The effect is that if the three terms 
in A4aa interfered maximally throughout phase space (which they do not), we would have 

(^{gg gggg; s) ■ a{gg gg-/-/; s) : a{gg ^ 7777; s) = 80 : 16 : 3 (35) 

while with no interference the numbers above would be 64 : 16 : 1. Thus the ratio of 



(7{gg 1111) to the total in Eq. (34) is 1/81 without interference, while if interference is 
maximal everywhere in phase space, the ratio is 1/33. In most regions of phase space, one 
of the three terms in M.aa will dominate, so interference effects will be small. But to be 
maximally conservative, since we have not performed the computation, we take the ratio 
1/33 for our upper bound. A full computation (or even a more detailed argument using the 
power-law dependence of F) would probably lead to a bound a factor of 1.5 to 2 stronger. 
This gives bounds on pp 7777 which are at least 33 times stronger than obtained just 



from Eq. (22), reducing the allowable 4-photon cross-sections to less than 25 femtobarns, as 



shown in Table III In particular, the case of A near 2, where the bound in [TO] was weakest, 
is where the unitarity bound is the strongest, below 0.15 fb. 

As we noted, this is obtained through a very conservative method, assuming (contrary 
to fact) that interference is maximal everywhere. Moreover, the reduction factor of 33 is 
increased to something closer to 40 by QCD corrections and by including processes involving 
Z bosons, such as ggZZ,gg'jZ, etc., in the final states. It would grow further if O also 
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couples to SU{2) gauge bosons. For these reasons we view 10 fb as a more likely bound. 
It is also worth noting that, were the bound saturated, requiring A ~ 1.2 and Ai = A™" 
as given in Table [11} then jet-plus-MET signals would significantly exceed Standard Model 
backgrounds at the LHC, giving a possible alternative discovery channel. 



B. General bounds 



The above situation is fairly generic. There is no reason to expect that any one pro- 
cess, especially one as experimentally attractive as gg — > 7777, dominates over all others. 
However, different processes cannot generically be combined together without additional 
calculation. For example, if Oi 7^ O2 and Ai 7^ A2, as we considered in most of this paper, 
then the choice of Smax for gg gggg is not the same as for gg 7777, and so their 
bounds are not simply related. Furthermore, although the four-gluon process is enhanced 
by color factors, it is proportional to a different three-point coefficient; C122 might be larger 
than Cm, and the indeed the latter could even be zero. In fact, we have implicitly assumed 
Cm = in our main discussion, because a non-zero value would give a stronger bound. 

The strongest model-independent bound we can obtain — using the unitarity constraints 
we have discussed above — is one given by assuming that the only large process at the scale 
Smax is gg 7777- This is in principle possible when Oi 7^ O2, so that A2 7^ Ai and 
Ai 7^ A2 in general. 

Our bounds in this more general setting, for various choices of Ai and A2, are shown in 



Table |IV| Interestingly, because our bounds on Ai are strong at low Ai but Smax is largest 
at higher Ai + 2A2, the bounds do not vary as widely as a function of Aj as one might have 
imagined. Note that for those values of Ai,A2 where the conformality constraint is more 
important than the experimental bound from jet-plus-MET, our bound depends only on 
Ai + 2A2; although Ai depends on Ai and A2 separately, the conformality constraint and 
the total cross-section a{gg {X}) both depend on A^^\ so that this dependence cancels 
out of our limit. 

Our bounds are smooth as the Aj pass through 2. This is because only the imaginary 
part of the unparticle two-point functions arises in our caculations. As a result, none of our 
intermediate steps require renormalization at A^ = 2. Conversely, note that we have cut off 
our table at A2 = 1.05. Although our bound would formally become still weaker as A2 — 1, 
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there is a separate constraint in this region. For A2 = 1, O2 is a free field [IHj, satisfying the 
Klein-Gordon equation, and therefore the OPE coefficient C122 ^ as A2 ^ 1 (with the 
unique exception of the case where Oi = {02)^, but then C122 1 and Ai ^ 2 so the rate 
cannot be large.) Consequently the four-photon production cross-section generated through 
the three-point function {O1O2O2) must be small as A2 1. 

Even though we are considering a much larger class of models, the limits we obtain are 
much stronger than those quoted in [10], especially at high Ai, A2. (For Oi = O2, as in [TO] , 
but generalizing by allowing Ai 7^ A2, the constraints are given along the diagonal, and are 
always below 1 pb.) However, we note that our bounds for Ai ~ 1.5, A2 ~ 1 - were they 
saturated - would still represent cross-sections of considerable phenomenological interest. 
One might have up to a few hundred events in the first year of running at the LHC. 

It is worth noting that where the bounds for pp 1111 lie well below 100 fb or so, 
this channel might not be the discovery channel. For the values of Ai shown in Table [IT], 
and for Ai < 1.4, the rate for jet-plus-MET at the LHC (for jet pt cuts of 250 GeV) is 



A2 

Ai 


1.05 


1.15 


1.25 


1.35 


1.45 


1.55 


1.65 


1.75 


1.85 


1.95 


1.05 


360 


170 


86 


45 


24 


13 


8 


5 


3 


2 


1.15 


1270 


640 


330 


180 


100 


58 


34 


21 


13 


8 


1.25 


2530 


1320 


720 


400 


230 


138 


83 


49 


27 


15 


1.35 


4270 


2330 


1120 


520 


250 


130 


66 


36 


20 


12 


1.45 


4020 


1690 


760 


360 


180 


91 


49 


27 


15 


9 


1.55 


2580 


1120 


520 


250 


126 


66 


36 


20 


12 


7 


1.65 


1690 


760 


360 


180 


91 


49 


27 
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TABLE IV: The maximum allowed values, in femtobarns, of the cross-section for pp 1111, as a 



function of Ai and A2, assuming Oi and O2 are different operators. (See Table III for the stronger 
bounds that apply \i Oi = 02-) Note that when the condition on Ai comes from the constraint of 
conformality, the bound depends only Ai + 2A2. 
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generally in the few pb range. This is somewhat larger than the standard model rate. Even 
though this measurement will be challenging in the early days of a new hadron collider, with 
substantial systematic errors, such a large excess in this channel might be convincing. This 
means that discovery of the new sector may well occur through the jet-plus-MET channel. 
In particular, this would almost certainly be the case in the model of [TU], given the tight 



(yet conservative) bounds in Table III For larger Ai or larger Ai the excess in jet-plus-MET 
may not be measurable, but also the four-photon rate would be even further reduced. 

Before concluding, we should re-emphasize the logic of our argument. Our claim is that 
if the cross-section for this process exceeds our bound, then conformal invariance must be 
strongly violated, which means that the universality of the "unparticle" dynamics is lost, 
and the calculations of [lOj, which assumed conformal invariance, are not valid. Instead, the 
production rate, and the kinematic distribution, would become highly model dependent. 

But we should hasten to add that large four-photon rates from a more general hidden 



sector are still possible. The bounds in Table [IV] only constrain a conformally invariant 
hidden sector. A large four-photon signal could come from other, non- "unparticle" hidden 
sectors — in particular from hidden valleys, which might or might not be conformal at 
high energy, but at low energy have strongly-broken conformal invariance and a mass gap. 
Examples of such theories are given in [121 EH] • Consequently, the four-photon experimental 
search channel, along with other multi-particle search channels, is of considerable interest 
in any case, and should be pursued model-independently. However, kinematic distributions 
will be very different from those in (TOj [161 (13) ^ind are highly model-dependent. 



C. Uncertainties on the bounds 



Our bounds, as they are upper bounds, do not need to account for any experimental 
considerations, such as triggering rates, acceptance or efficiency, event selection cuts and 
the like, which can only reduce the number of events. Indeed such considerations enter only 
in our determination of A™" from existing experimental data. Because the gg 7777 
cross-section is largest at large s, giving four photons which typically have momenta in the 
few hundred GeV range, neither triggering, efficiency or even geometric acceptance are likely 
to reduce significantly the number of observed events at the LHC This is especially true if 
a loose criterion (such as demanding only three of the four photons be observed) is applied 
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in the analysis. 

Still, our results have multiple sources of uncertainties. For example, we ignored K-factors 
which would have given us a stronger bound on Ai, but which also would have given us a 
larger cross-section for gg —>■ {X} and therefore a weaker bound on pp — > 7777; these effects 
most likely cancel to a good approximation. We also did not use the most updated parton 
distribution functions, and in any case applied them only in a leading order approximation. 
We neglected some experimental efficiencies in our extraction of Ai, but were conservative 
in our use of the CDF data from |25]. We included only the largest jet-plus-MET process 
at the Tevatron, worked only at leading order, and treated errors in the CDF data using 
crude estimates of systematic and statistical errors. Also we have used results from only 1.1 
inverse fb; unpublished limits have probably improved somewhat. 

But the dominant source of uncertainty in our bound comes from our choice of the 
parameter ( defining Smax, and for this reason it does not make sense for us to reduce the 



uncertainties mentioned in the previous paragraph. We chose to use C = 3 in Eq- (17). 
Using C = ^ could loosen our bounds by a factor of order 3-5. On the other hand, such a 
choice puts the peak cross-section right at the value of s where the unitarity bound is kicking 
in, which means that conformal invariance is breaking down precisely where a prediction is 
most needed. One could also argue that C = | is a better choice, which would tighten the 
bounds by a factor of order 2. In any case, one must view this choice as one of taste. But in 
addition we think it highly unlikely that a strict unitarity bound would be fully saturated 
in any physical model. It is much more probable that either conformal invariance will break 
down below Smax, or that the pointlike interaction between the two sectors will develop a 
form factor below Smax- Thus we expect that typically a breakdown of the methods of 110] 
occurs well below the energy where the gg —>■ 7777 cross-section formally would exceed the 
gg {X} cross-section. In this sense, we expect that our bounds, though imprecise, are 
actually quite conservative. 



VII. COMMENTS ON OTHER MULTI-PARTICLE PROCESSES 

There are many other processes to which this type of unitarity bound should be applied, 
each with its own features which we did not fully explore here. In particular, this type of 
bound is powerful whenever the couplings between the two sectors are non-renormalizable. 
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a condition which ensures that a process such as gg — >■ 7777 grows with energy relative to 
gg —>■ {X}. (Actually it is enough that the couplings involving the final state particles, in 
our case 02F^yF^^ , be non-renormalizable.) 

An example where our bound would not be strong is in the process qq —>■ 
through three scalar operators of A ~ 1, as considered in [TB]. Here the operator coupling 
the two sectors (after the Higgs gets an expectation value) has dimension near four if the Aj 
are not far above 1. But conversely, as was demonstrated in [16], the lack of rapid growth at 
high energy also means there is no suppression at low energy, and therefore Tevatron limits 
are very strong. Meanwhile, our arguments do apply if the Aj are significantly larger than 
1. 



We argued in Sec. IIIB, however, that this case is not physically reasonable anyway. 
Large fiavor-changing neutral currents are essentially impossible to avoid if one couples a 
new sector through chirality-fiipping operators (as would be the case for scalars) to light 
quarks and leptons. 

The problem of fiavor-changing currents would be alleviated in models where the cou- 
plings to the quarks and leptons are weighted by mass, so that no additional fiavor dynamics 
is introduced. In this case one might consider gg bbbb or gg t^t^t^t^ . Here the 
bounds from our methods would be weak. Fermilab production of this process would not be 
strongly constrained in the case of bbbb. However the trilepton searches at Fermilab would 
significantly constrain the four-tau final state. Another possibility would involve gg — > 7766 
or gg 'y-ftt. Our bound for the sum of these processes is roughly 30 times weaker than for 
gg —* '~f'~f'~f'~f. Backgrounds of course are larger too, but limits from Fermilab on 7766 may 
be rather weak, and on 77^ will be very limited because of kinematic constraints and low 
statistics. This case might merit additional exploration. 

Another possibility involves couplings of standard model particles to non-scalar operators 
in the conformal field theory. In some cases the couplings to light quarks and leptons would 
be chirality preserving and need not introduce any new fiavor dependence. Because unitarity 
requires vector operators have dimension 3 or greater, and tensor operators to have dimension 
4 or greater, their couplings to the standard model are always non-renormalizable. Four- 
particle final states generated through vector operators have growing cross-sections. This 
means Tevatron bounds on processes such as qq —>■ via vectors operators are weak, 

but conversely our unitarity constraints are very strong. 
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For example, one option with no fine-tuning might involve the possibility of a three-point 
function involving two pseudoscalar operators and a vector operator. Consider the process 
99 77^"*"^" which would arise in a theory which has, in addition to the two couplings in 



Eq. (15), a third coupling 

-^V, E^a^E, (36) 

2 i 

where Ei is a left-handed charged antilepton e+,/i"'",r+. Because the vector operator 
must have dimension Ay > 3, the constraints obtained via our methods are 10-30 times 
stronger than those for gg — >■ 7777, with the maximum allowed cross-section being of order 
100 fb. 



VIII. CONCLUSIONS AND OUTLOOK 



We considered an example of a multi-particle process mediated by a hidden sector that 
is conformally invariant, along the lines of [in]. Conformal invariance makes the process 
predictable, in a way that depends only on the dimensions of the operators, up to an 
overall normalization. We have shown that the total cross-sections for such processes are 
strongly constrained by requiring both conformal invariance and unitarity. The constraint 
is generally stronger when the products of standard-model and hidden-sector operators that 
appear in the action have dimensions significantly larger than 4. This is because such non- 
renormalizable interactions generate cross-sections that grow rapidly with energy, and will 
become larger than the total hidden-sector production cross-section at an energy that is of 
order A, the scale of the coupling of the two sectors. 

In particular, we saw that, in the model suggested by [10], the process gg 7777 is 
constrained to lie below 25 fb. Moreover, for operators with dimension A < 1.5, saturating 
this bound would require a scale A so low that the rate for jet-plus-MET would be larger, 
even at moderate pt-, than the standard model rate. For operators with A > 1.5, the bound 
on gg 7777 is below 3 fb. 

However, relaxing the restrictive conditions in [lOj allowed us to raise the limits on the 
four-photon cross-section, giving substantial LHC signals potentially as large as a few pb. 
But we emphasize that we believe that this is only the beginning of the story. More sophis- 
ticated constraints from unitarity appear possible. If so, the quantitative results obtained 
here will be tightened further. We hope to report on this, and clarify the phenomenological 
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situation, in a subsequent publication. 

As we noted, our methods apply more widely. Processes such as gg —>■ '-f-fbb with scalar 
operators coupling to heavy flavor fermion-bilinears, which grow more slowly with energy 
than gg 7777, may be less constrained by unitarity, while processes involving vector- 
operators, such as gg 'j'yi^i' , which grow more rapidly, are more constrained. However, 
experimental constraints from Fermilab are stronger in the former case than the latter, 
precisely because of this difference in energy dependence. 

Our quantitative results do suffer from some ambiguities. On the one hand, we have been 
very conservative in our numbers. We believe that realistic limits are at least a factor of 2 
or 3 stronger than we have claimed. Also, in real models the bounds that we obtained will 
rarely be saturated, and even when they are, it is unlikely that the process which saturates 
the bound will be the easiest to observe, as gg —>■ 7777 would be. On the other hand, one 
could take an even more conservative view regarding our definition of Smax, and get bounds 
weaker by a factor of 3 or so. However there is no way to weaken our bounds by much more 
than this, except by giving up conformal invariance, and with it the model-independent 
predictions of the "unparticle" scenario. 
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